Next-generation sequencing reveals mitogenome diversity in plasma extracellular vesicles from colorectal cancer patients.

BACKGROUND: Recent reports have demonstrated that the entire mitochondrial genome can be secreted in extracellular vesicles (EVs), but the biological attributes of this cell-free mitochondrial DNA (mtDNA) remain insufficiently understood. We used next-generation sequencing to compare plasma EV-derived mtDNA to that of whole blood (WB), peripheral blood mononuclear cells (PBMCs), and formalin-fixed paraffin-embedded (FFPE) tumor tissue from eight rectal cancer patients and WB and fresh-frozen (FF) tumor tissue from eight colon cancer patients. METHODS: Total DNA was isolated before the mtDNA was enriched by PCR with either two primer sets generating two long products or multiple primer sets (for the FFPE tumors), prior to the sequencing. mtDNA diversity was assessed as the total variant number, level of heteroplasmy (mutant mtDNA copies mixed with wild-type copies), variant distribution within the protein-coding genes, and the predicted functional effect of the variants in the different sample types. Differences between groups were compared by paired Student's t-test or ANOVA with Dunnett's multiple comparison tests when comparing matched samples from patients. Mann-Whitney U test was used when comparing differences between the cancer types and patient groups. Pearson correlation analysis was performed. RESULTS: In both cancer types, EV mtDNA presented twice as many variants and had significantly more low-level heteroplasmy than WB mtDNA. The EV mtDNA variants were clustered in the coding regions, and the proportion of EV mtDNA variants that were missense mutations (i.e., estimated to moderately affect the mitochondrial protein function) was significantly higher than in WB and tumor tissues. Nonsense mutations (i.e., estimated to highly affect the mitochondrial protein function) were only observed in the tumor tissues and EVs. CONCLUSION: Taken together, plasma EV mtDNA in CRC patients exhibits a high degree of diversity. TRIAL REGISTRATION: ClinicalTrials.gov: NCT01816607 . Registered 22 March 2013.

Longitudinal whole blood transcriptomic analysis characterizes neutrophil activation and interferon signaling in moderate and severe COVID-19.

A maladaptive inflammatory response has been implicated in the pathogenesis of severe COVID-19. This study aimed to characterize the temporal dynamics of this response and investigate whether severe disease is associated with distinct gene expression patterns. We performed microarray analysis of serial whole blood RNA samples from 17 patients with severe COVID-19, 15 patients with moderate disease and 11 healthy controls. All study subjects were unvaccinated. We assessed whole blood gene expression patterns by differential gene expression analysis, gene set enrichment, two clustering methods and estimated relative leukocyte abundance using CIBERSORT. Neutrophils, platelets, cytokine signaling, and the coagulation system were activated in COVID-19, and this broad immune activation was more pronounced in severe vs. moderate disease. We observed two different trajectories of neutrophil-associated genes, indicating the emergence of a more immature neutrophil phenotype over time. Interferon-associated genes were strongly enriched in early COVID-19 before falling markedly, with modest severity-associated differences in trajectory. In conclusion, COVID-19 necessitating hospitalization is associated with a broad inflammatory response, which is more pronounced in severe disease. Our data suggest a progressively more immature circulating neutrophil phenotype over time. Interferon signaling is enriched in COVID-19 but does not seem to drive severe disease.

TiSA: TimeSeriesAnalysis-a pipeline for the analysis of longitudinal transcriptomics data.

Improved transcriptomic sequencing technologies now make it possible to perform longitudinal experiments, thus generating a large amount of data. Currently, there are no dedicated or comprehensive methods for the analysis of these experiments. In this article, we describe our TimeSeries Analysis pipeline (TiSA) which combines differential gene expression, clustering based on recursive thresholding, and a functional enrichment analysis. Differential gene expression is performed for both the temporal and conditional axes. Clustering is performed on the identified differentially expressed genes, with each cluster being evaluated using a functional enrichment analysis. We show that TiSA can be used to analyse longitudinal transcriptomic data from both microarrays and RNA-seq, as well as small, large, and/or datasets with missing data points. The tested datasets ranged in complexity, some originating from cell lines while another was from a longitudinal experiment of severity in COVID-19 patients. We have also included custom figures to aid with the biological interpretation of the data, these plots include Principal Component Analyses, Multi Dimensional Scaling plots, functional enrichment dotplots, trajectory plots, and complex heatmaps showing the broad overview of results. To date, TiSA is the first pipeline to provide an easy solution to the analysis of longitudinal transcriptomics experiments.

Krill oil protects dopaminergic neurons from age-related degeneration through temporal transcriptome rewiring and suppression of several hallmarks of aging.

There is accumulating evidence that interfering with the basic aging mechanisms can enhance healthy longevity. The interventional/therapeutic strategies targeting multiple aging hallmarks could be more effective than targeting one hallmark. While health-promoting qualities of marine oils have been extensively studied, the underlying molecular mechanisms are not fully understood. Lipid extracts from Antarctic krill are rich in long-chain omega-3 fatty acids choline, and astaxanthin. Here, we used C. elegans and human cells to investigate whether krill oil promotes healthy aging. In a C. elegans model of Parkinson´s disease, we show that krill oil protects dopaminergic neurons from aging-related degeneration, decreases alpha-synuclein aggregation, and improves dopamine-dependent behavior and cognition. Krill oil rewires distinct gene expression programs that contribute to attenuating several aging hallmarks, including oxidative stress, proteotoxic stress, senescence, genomic instability, and mitochondrial dysfunction. Mechanistically, krill oil increases neuronal resilience through temporal transcriptome rewiring to promote anti-oxidative stress and anti-inflammation via healthspan regulating transcription factors such as SNK-1. Moreover, krill oil promotes dopaminergic neuron survival through regulation of synaptic transmission and neuronal functions via PBO-2 and RIM-1. Collectively, krill oil rewires global gene expression programs and promotes healthy aging via abrogating multiple aging hallmarks, suggesting directions for further pre-clinical and clinical explorations.

Mucosal Gene Transcript Signatures in Treatment Naïve Inflammatory Bowel Disease: A Comparative Analysis of Disease to Symptomatic and Healthy Controls in the European IBD-Character Cohort.

BACKGROUND: Studies of the mucosal transcriptomic landscape have given new insight into the pathogenesis of inflammatory bowel disease (IBD). Recently, the predictive biomarker potential of gene expression signatures has been explored. To further investigate the mucosal gene expression in IBD, we recruited a cohort of treatment naïve patients and compared them to both symptomatic and healthy controls. METHODS: Altogether, 323 subjects were included: Crohn's disease (N = 75), ulcerative colitis (N = 87) and IBD unclassified (N = 3). Additionally, there were two control groups: symptomatic controls (N = 131) and healthy controls (N = 27). Mucosal biopsies were collected during ileocolonoscopy and gene expression in inflamed and non-inflamed mucosa was explored. Gene expression profiling was performed using Agilent G3 Human Gene Expression 860K v3 One-Color microarray. We recorded information about treatment escalation to anti-TNF agents or surgery, and anti-TNF response, to explore predictive opportunities of the mucosal transcriptome. RESULTS: Gene expression profiles in symptomatic controls in whom IBD had been excluded resembled that of IBD patients and diverged from that of healthy controls. In non-inflamed Crohn's disease and ulcerative colitis, gene set enrichment analysis revealed dysregulation of pathways involved in basic cellular biological processes. Mitochondria-associated pathways were dysregulated both in non-inflamed and inflamed Crohn's disease and ulcerative colitis (>2.6 normalized enrichment scores <-1.8). Gene expression signatures of Crohn's disease and ulcerative colitis did not predict time for treatment escalation (p = 0.175). No significant association was found between gene expression signatures and anti-TNF response. CONCLUSION: Non-inflamed samples are probably superior to inflamed samples when exploring gene expression signatures in IBD and might reveal underlying mechanisms central for disease initiation. The gene expression signatures of the control groups were related to if they were symptomatic or not, which may have important implications for future study designs.

ZBTB11 dysfunction: spectrum of brain abnormalities, biochemical signature and cellular consequences.

Bi-allelic pathogenic variants in ZBTB11 have been associated with intellectual developmental disorder, autosomal recessive 69 (MRT69; OMIM 618383). We report five patients from three families with novel, bi-allelic variants in ZBTB11. We have expanded the clinical phenotype of MRT69, documenting varied severity of atrophy affecting different brain regions and described combined malonic and methylmalonic aciduria as a biochemical manifestation. As ZBTB11 encodes for a transcriptional regulator, we performeded chromatin immunoprecipitation-sequencing targeting ZBTB11 in fibroblasts from patients and controls. Chromatin immunoprecipitation-sequencing revealed binding of wild-type ZBTB11 to promoters in 238 genes, among which genes encoding proteins involved in mitochondrial functions and RNA processing are over-represented. Mutated ZBTB11 showed reduced binding to 61 of the targeted genes, indicating that the variants act as loss of function. Most of these genes are related to mitochondrial functions. Transcriptome analysis of the patient fibroblasts revealed dysregulation of mitochondrial functions. In addition, we uncovered that reduced binding of the mutated ZBTB11 to ACSF3 leads to decreased ACSF3 transcript level, explaining combined malonic and methylmalonic aciduria. Collectively, these results expand the clinical spectrum of ZBTB11-related neurological disease and give insight into the pathophysiology in which the dysfunctional ZBTB11 affect mitochondrial functions and RNA processing contributing to the neurological and biochemical phenotypes.

Comparable cancer-relevant mutation profiles in synchronous ductal carcinoma in situ and invasive breast cancer.

BACKGROUND: Ductal carcinoma in situ (DCIS) comprises a diverse group of preinvasive lesions in the breast and poses a considerable clinical challenge due to lack of markers of progression. Genomic alterations are to a large extent similar in DCIS and invasive carcinomas, although differences in copy number aberrations, gene expression patterns, and mutations exist. In mixed tumors with synchronous invasive breast cancer (IBC) and DCIS, it is still unclear to what extent invasive tumor cells are directly derived from the DCIS cells. AIM: Our aim was to compare cancer-relevant mutation profiles of different cellular compartments in mixed DCIS/IBC and pure DCIS tumors. METHODS AND RESULTS: We performed targeted sequencing of 50 oncogenes in microdissected tissue from three different epithelial cell compartments (in situ, invasive, and normal adjacent epithelium) from 26 mixed breast carcinomas. In total, 44 tissue samples (19 invasive, 16 in situ, 9 normal) were subjected to sequencing using the Ion Torrent platform and the AmpliSeq Cancer Hotspot Panel v2. For comparison, 10 additional, pure DCIS lesions were sequenced. Across all mixed samples, we detected 23 variants previously described in cancer. The most commonly affected genes were TP53, PIK3CA, and ERBB2. The PIK3CA:p.H1047R variant was found in nine samples from six patients. Most variants detected in invasive compartments were also found in the corresponding in situ cell compartment indicating a clonal relationship between the tumor stages. A lower frequency of variants were observed in pure DCIS lesions. CONCLUSION: Similar mutation profiles between in situ and invasive cell compartments indicate a similar origin of the two tumor stages in mixed breast tumors. The lower number of potential driver variants found in pure DCIS compared with the in situ cell compartments of mixed tumors may imply that pure DCIS is captured earlier in the path of progression to invasive disease.

Loss of progesterone receptor is associated with distinct tyrosine kinase profiles in breast cancer.

PURPOSE: The aim of this study was to assess protein tyrosine kinase profiles in primary breast cancer samples in correlation with the distinct hormone and growth receptor profiles ER, PR, and HER2. EXPERIMENTAL DESIGN: Pamchip® microarrays were used to measure the phosphorylation of 144 tyrosine kinase substrates in 29 ER+ breast cancer samples and cell lines MCF7, BT474 and ZR75-1. mRNA expression data from the METABRIC cohort and publicly available PR chip-sequencing data were used for validation purposes, together with RT-PCR. RESULTS: In ER+ breast tumors and cell lines, we observed that the loss of PR expression correlated to higher kinase activity in samples and cell lines that were HER2-. A number of kinases, representing mostly proteins within the PI3K/AKT pathway, were identified as responsible for the differential phosphorylation between PR- and PR+ in ER+/HER2- tumors. We used the METABRIC cohort to analyze mRNA expression from 977 ER+/HER2- breast cancers. Twenty four kinase-encoding genes were identified as differentially expressed between PR+ and PR-, dividing ER+/HER2- samples in two distinct clusters with significant differences in survival (p < 0.05). Four kinase genes, LCK, FRK, FGFR4, and MST1R, were identified as potential direct targets of PR. CONCLUSIONS: Our results suggest that the PR status has a profound effect on tyrosine kinases, especially for FGFR4 and LCK genes, in ER+/HER2- breast cancer patients. The influence of these genes on the PI3K/AKT signaling pathway may potentially lead to novel drug targets for ER+/PR- breast cancer patients.

The NEOLETEXE trial: a neoadjuvant cross-over study exploring the lack of cross resistance between aromatase inhibitors.

The aromatase inhibitor letrozole (Femar®/Femara®) and the aromatase inactivator exemestane (Aromasin®) differ in their biochemical effect on the aromatase enzyme. Letrozole is a competitive aromatase inhibitor while exemestane binds irreversibly to the aromatase enzyme. This pharmacological difference is of clinical interest since a lack of cross-resistance has been documented. It has been demonstrated in several clinical trials that exemestane may cause a disease regression following resistance to nonsteroidal aromatase inhibitors. The exact mechanism(s) behind this phenomenon is yet unknown. Here, we present the NEOLETEXE trial with the aim of exploring the individual mechanisms involved behind the observed lack of cross resistance. Clinical trial registration: The trial has been approved by the Regional Ethics Committee of South-East Norway (project number 2015/84).

Breast cancer quantitative proteome and proteogenomic landscape.

In the preceding decades, molecular characterization has revolutionized breast cancer (BC) research and therapeutic approaches. Presented herein, an unbiased analysis of breast tumor proteomes, inclusive of 9995 proteins quantified across all tumors, for the first time recapitulates BC subtypes. Additionally, poor-prognosis basal-like and luminal B tumors are further subdivided by immune component infiltration, suggesting the current classification is incomplete. Proteome-based networks distinguish functional protein modules for breast tumor groups, with co-expression of EGFR and MET marking ductal carcinoma in situ regions of normal-like tumors and lending to a more accurate classification of this poorly defined subtype. Genes included within prognostic mRNA panels have significantly higher than average mRNA-protein correlations, and gene copy number alterations are dampened at the protein-level; underscoring the value of proteome quantification for prognostication and phenotypic classification. Furthermore, protein products mapping to non-coding genomic regions are identified; highlighting a potential new class of tumor-specific immunotherapeutic targets.

Noninvasive profiling of serum cytokines in breast cancer patients and clinicopathological characteristics.

Cancers elicit an immune response by modifying the microenvironment. The immune system plays a pivotal role in cancer recognition and eradication. While the potential clinical value of infiltrating lymphocytes at the tumor site has been assessed in breast cancer, circulating cytokines - the molecules coordinating and fine-tuning immune response - are still poorly characterized. Using two breast cancer cohorts (MicMa, n = 131, DCTB, n = 28) and the multiplex Luminex platform, we measured the levels of 27 cytokines in the serum of breast cancer patients prior to treatment. We investigated the cytokine levels in relation to clinicopathological characteristics and in perspective of the tumor infiltrating immune cells predicted from the bulk mRNA expression data. Unsupervised clustering analysis of the serum cytokine levels in the MicMa cohort identified a cluster of pro-inflammatory, pro-angiogenic, and Th2-related cytokines which was associated with poor prognosis. Notably high levels of platelet derived growth factor BB (PDGF) reflected a more aggressive tumor phenotype and larger tumor size. A significant positive correlation between serum levels of interferon gamma-induced protein 10 (IP10) and its mRNA expression at the tumor site suggested that tumor-IP10-production may outflow to the bloodstream. High IP10 serum levels were associated with a worse prognosis. Finally, we found serum levels of both PDGF and IP10 associated with enrichment scores of specific tumor infiltrating immune cells. Our study suggests that monitoring cytokine circulating levels in breast cancer could be used to characterize breast cancers and the immune composition of their microenvironment through readily available biological material.

Serum cytokine levels in breast cancer patients during neoadjuvant treatment with bevacizumab.

A high concentration of circulating vascular endothelial growth factor (VEGF) in cancer patients is associated with an aggressive tumor phenotype. Here, serum levels of 27 cytokines and blood cell counts were assessed in breast cancer patients receiving neoadjuvant chemotherapy with or without bevacizumab (Bev) in a randomized cohort of 132 patients with non-metastatic HER2-negative tumors. Cytokine levels were determined prior to treatment and at various time-points. The cytotoxic chemotherapy regimen of fluorouracil, epirubicin, and cyclophosphamide (FEC) had a profound impact on both circulating white blood cells and circulating cytokine levels. At the end of FEC treatment, the global decrease in cytokine levels correlated with the drop in white blood cell counts and was significantly greater in the patients of the Bev arm for cytokines, such as VEGF-A, IL-12, IP-10 and IL-10. Among patients who received Bev, those with pathological complete response (pCR) exhibited significantly lower levels of VEGF-A, IFN-γ, TNF-α and IL-4 than patients without pCR. This effect was not observed in the chemotherapy-only arm. Certain circulating cytokine profiles were found to correlate with different immune cell types at the tumor site. For the Bev arm patients, the serum cytokine levels correlated with higher levels of cytotoxic T cells at the end of the therapy regimen, which was indicative of treatment response. The higher response rate for Bev-treated patients and stronger correlations between serum cytokine levels and infiltrating CD8T cells merits further investigation.

Age, estrogen, and immune response in breast adenocarcinoma and adjacent normal tissue.

Chronic inflammation promotes breast tumor growth and invasion by accelerating angiogenesis and tissue remodeling in the tumor microenvironment. There is a complex relationship between inflammation and estrogen, which drives the growth of 70 percent of breast tumors. While low levels of estrogen exposure stimulate macrophages and other inflammatory cell populations, very high levels are immune suppressive. Breast tumor incidence is increased by obesity and age, which interact to influence inflammatory cell populations in normal breast tissue. To characterize the impact of these factors on tumors and the tumor microenvironment, we measured gene expression in 195 breast adenocarcinomas and matched adjacent normal breast tissue samples collected at Akershus University Hospital (AHUS). Age and Body Mass Index (BMI) were independently associated with inflammation in adjacent normal tissue but not tumors. Estrogen Receptor (ER)-negative tumors had elevated macrophage expression compared with matched normal tissue, but ER-positive tumors showed an unexpected decrease in macrophage expression. We found an inverse relationship between the increase in tumor estrogen pathway expression compared with adjacent normal tissue and tumor macrophage score. We validated this finding in 126 breast tumor-normal pairs from the previously published METABRIC cohort. We developed a novel statistic, the Rewiring Coefficient, to quantify the rewiring of gene co-expression networks at the level of individual genes. Differential correlation analysis demonstrated distinct pathways were rewired during tumorigenesis. Our data support an immune suppressive effect of high doses of estrogen signaling in breast tumor microenvironment, suggesting that this effect contributes to the greater presence of prognostic and therapeutically relevant immune cells in ER-negative tumors.

Integrative clustering reveals a novel split in the luminal A subtype of breast cancer with impact on outcome.

BACKGROUND: Breast cancer is a heterogeneous disease at the clinical and molecular level. In this study we integrate classifications extracted from five different molecular levels in order to identify integrated subtypes. METHODS: Tumor tissue from 425 patients with primary breast cancer from the Oslo2 study was cut and blended, and divided into fractions for DNA, RNA and protein isolation and metabolomics, allowing the acquisition of representative and comparable molecular data. Patients were stratified into groups based on their tumor characteristics from five different molecular levels, using various clustering methods. Finally, all previously identified and newly determined subgroups were combined in a multilevel classification using a "cluster-of-clusters" approach with consensus clustering. RESULTS: Based on DNA copy number data, tumors were categorized into three groups according to the complex arm aberration index. mRNA expression profiles divided tumors into five molecular subgroups according to PAM50 subtyping, and clustering based on microRNA expression revealed four subgroups. Reverse-phase protein array data divided tumors into five subgroups. Hierarchical clustering of tumor metabolic profiles revealed three clusters. Combining DNA copy number and mRNA expression classified tumors into seven clusters based on pathway activity levels, and tumors were classified into ten subtypes using integrative clustering. The final consensus clustering that incorporated all aforementioned subtypes revealed six major groups. Five corresponded well with the mRNA subtypes, while a sixth group resulted from a split of the luminal A subtype; these tumors belonged to distinct microRNA clusters. Gain-of-function studies using MCF-7 cells showed that microRNAs differentially expressed between the luminal A clusters were important for cancer cell survival. These microRNAs were used to validate the split in luminal A tumors in four independent breast cancer cohorts. In two cohorts the microRNAs divided tumors into subgroups with significantly different outcomes, and in another a trend was observed. CONCLUSIONS: The six integrated subtypes identified confirm the heterogeneity of breast cancer and show that finer subdivisions of subtypes are evident. Increasing knowledge of the heterogeneity of the luminal A subtype may add pivotal information to guide therapeutic choices, evidently bringing us closer to improved treatment for this largest subgroup of breast cancer.

Subtype-specific micro-RNA expression signatures in breast cancer progression.

Robust markers of invasiveness may help reduce the overtreatment of in situ carcinomas. Breast cancer is a heterogeneous disease and biological mechanisms for carcinogenesis vary between subtypes. Stratification by subtype is therefore necessary to identify relevant and robust signatures of invasive disease. We have identified microRNA (miRNA) alterations during breast cancer progression in two separate datasets and used stratification and external validation to strengthen the findings. We analyzed two separate datasets (METABRIC and AHUS) consisting of a total of 186 normal breast tissue samples, 18 ductal carcinoma in situ (DCIS) and 1,338 invasive breast carcinomas. Validation in a separate dataset and stratification by molecular subtypes based on immunohistochemistry, PAM50 and integrated cluster classifications were performed. We propose subtype-specific miRNA signatures of invasive carcinoma and a validated signature of DCIS. miRNAs included in the invasive signatures include downregulation of miR-139-5p in aggressive subtypes and upregulation of miR-29c-5p expression in the luminal subtypes. No miRNAs were differentially expressed in the transition from DCIS to invasive carcinomas on the whole, indicating the need for subtype stratification. A total of 27 miRNAs were included in our proposed DCIS signature. Significant alterations of expression included upregulation of miR-21-5p and the miR-200 family and downregulation of let-7 family members in DCIS samples. The signatures proposed here can form the basis for studies exploring DCIS samples with increased invasive potential and serum biomarkers for in situ and invasive breast cancer.

Gene expression analysis supports tumor threshold over 2.0 cm for T-category breast cancer.

Tumor size, as indicated by the T-category, is known as a strong prognostic indicator for breast cancer. It is common practice to distinguish the T1 and T2 groups at a tumor size of 2.0 cm. We investigated the 2.0-cm rule from a new point of view. Here, we try to find the optimal threshold based on the differences between the gene expression profiles of the T1 and T2 groups (as defined by the threshold). We developed a numerical algorithm to measure the overall differential gene expression between patients with smaller tumors and those with larger tumors among multiple expression datasets from different studies. We confirmed the performance of the proposed algorithm by a simulation study and then applied it to three different studies conducted at two Norwegian hospitals. We found that the maximum difference in gene expression is obtained at a threshold of 2.2-2.4 cm, and we confirmed that the optimum threshold was over 2.0 cm, as indicated by a validation study using five publicly available expression datasets. Furthermore, we observed a significant differentiation between the two threshold groups in terms of time to local recurrence for the Norwegian datasets. In addition, we performed an associated network and canonical pathway analyses for the genes differentially expressed between tumors below and above the given thresholds, 2.0 and 2.4 cm, using the Norwegian datasets. The associated network function illustrated a cellular assembly of the genes for the 2.0-cm threshold: an energy production for the 2.4-cm threshold and an enrichment in lipid metabolism based on the genes in the intersection for the 2.0- and 2.4-cm thresholds.

Canine Mammary Tumours Are Affected by Frequent Copy Number Aberrations, including Amplification of MYC and Loss of PTEN.

BACKGROUND: Copy number aberrations frequently occur during the development of many cancers. Such events affect dosage of involved genes and may cause further genomic instability and progression of cancer. In this survey, canine SNP microarrays were used to study 117 canine mammary tumours from 69 dogs. RESULTS: We found a high occurrence of copy number aberrations in canine mammary tumours, losses being more frequent than gains. Increased frequency of aberrations and loss of heterozygosity were positively correlated with increased malignancy in terms of histopathological diagnosis. One of the most highly recurrently amplified regions harbored the MYC gene. PTEN was located to a frequently lost region and also homozygously deleted in five tumours. Thus, deregulation of these genes due to copy number aberrations appears to be an important event in canine mammary tumour development. Other potential contributors to canine mammary tumour pathogenesis are COL9A3, INPP5A, CYP2E1 and RB1. The present study also shows that a more detailed analysis of chromosomal aberrations associated with histopathological parameters may aid in identifying specific genes associated with canine mammary tumour progression. CONCLUSIONS: The high frequency of copy number aberrations is a prominent feature of canine mammary tumours as seen in other canine and human cancers. Our findings share several features with corresponding studies in human breast tumours and strengthen the dog as a suitable model organism for this disease.

Glycan-related gene expression signatures in breast cancer subtypes; relation to survival.

Alterations in glycan structures are early signs of malignancy and have recently been proposed to be in part a driving force behind malignant transformation. Here, we explore whether differences in expression of genes related to the process of glycosylation exist between breast carcinoma subtypes - and look for their association to clinical parameters. Five expression datasets of 454 invasive breast carcinomas, 31 ductal carcinomas in situ (DCIS), and 79 non-malignant breast tissue samples were analysed. Results were validated in 1960 breast carcinomas. 419 genes encoding glycosylation-related proteins were selected. The DCIS samples appeared expression-wise similar to carcinomas, showing altered gene expression related to glycosaminoglycans (GAGs) and N-glycans when compared to non-malignant samples. In-situ lesions with different aggressiveness potentials demonstrated changes in glycosaminoglycan sulfation and adhesion proteins. Subtype-specific expression patterns revealed down-regulation of genes encoding glycan-binding proteins in the luminal A and B subtypes. Clustering basal-like samples using a consensus list of genes differentially expressed across discovery datasets produced two clusters with significantly differing prognosis in the validation dataset. Finally, our analyses suggest that glycolipids may play an important role in carcinogenesis of breast tumors - as demonstrated by association of B3GNT5 and UGCG genes to patient survival. In conclusion, most glycan-specific changes occur early in the carcinogenic process. We have identified glycan-related alterations specific to breast cancer subtypes including a prognostic signature for two basal-like subgroups. Future research in this area may potentially lead to markers for better prognostication and treatment stratification of breast cancer patients.

Long non-coding RNAs differentially expressed between normal versus primary breast tumor tissues disclose converse changes to breast cancer-related protein-coding genes.

Breast cancer, the second leading cause of cancer death in women, is a highly heterogeneous disease, characterized by distinct genomic and transcriptomic profiles. Transcriptome analyses prevalently assessed protein-coding genes; however, the majority of the mammalian genome is expressed in numerous non-coding transcripts. Emerging evidence supports that many of these non-coding RNAs are specifically expressed during development, tumorigenesis, and metastasis. The focus of this study was to investigate the expression features and molecular characteristics of long non-coding RNAs (lncRNAs) in breast cancer. We investigated 26 breast tumor and 5 normal tissue samples utilizing a custom expression microarray enclosing probes for mRNAs as well as novel and previously identified lncRNAs. We identified more than 19,000 unique regions significantly differentially expressed between normal versus breast tumor tissue, half of these regions were non-coding without any evidence for functional open reading frames or sequence similarity to known proteins. The identified non-coding regions were primarily located in introns (53%) or in the intergenic space (33%), frequently orientated in antisense-direction of protein-coding genes (14%), and commonly distributed at promoter-, transcription factor binding-, or enhancer-sites. Analyzing the most diverse mRNA breast cancer subtypes Basal-like versus Luminal A and B resulted in 3,025 significantly differentially expressed unique loci, including 682 (23%) for non-coding transcripts. A notable number of differentially expressed protein-coding genes displayed non-synonymous expression changes compared to their nearest differentially expressed lncRNA, including an antisense lncRNA strongly anticorrelated to the mRNA coding for histone deacetylase 3 (HDAC3), which was investigated in more detail. Previously identified chromatin-associated lncRNAs (CARs) were predominantly downregulated in breast tumor samples, including CARs located in the protein-coding genes for CALD1, FTX, and HNRNPH1. In conclusion, a number of differentially expressed lncRNAs have been identified with relation to cancer-related protein-coding genes.